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Quantum Fisher information is a key concept in the field of quantum metrology, which aims to 
enhance the accuracy of parameter estimation by using quantum resources. In this paper, utilizing 
a representation of quantum Fisher information for a general unitary parametrization process, we 
study unitary parametrization processes governed by su(2) dynamics. We obtain the analytical 
expression for the Hermitian operator of the parametrization and the maximal quantum Fisher 
information. We find that the maximal quantum Fisher information over the parameter space 
consists of two parts; one is quadratic in time and the other oscillates with time. We apply our 
result to the estimation of a magnetic field and obtained the maximal quantum Fisher information. 

We further discuss a driving field with a time-dependent Hamiltonian and find that the maximal 
quantum Fisher information of the driving frequency attains the optimum when it is in resonance 
with the atomic frequency. 

PACS numbers: 03.67.-a, 03.65.Ta, 06.20.-f 


I. INTRODUCTION 

The past two decades have witnessed a dramatic de¬ 
velopment of quantum metrology mug, which is rooted 
in the theory of quantum parameter estimation. The 
estimation of unknown parameters plays an important 
role in physics and engineering. By taking advantage of 
quantum resources such as entanglement and squeezing, 
quantum metrology promises a higher precision in pa¬ 
rameter estimation than its classical counterpart. There¬ 
fore many practical applications have been suggested to 
reap this benefit, including the detection of gravitational 
radiation hmd, quantum frequency standards IMD, 
and quantum imaging [22H24| . 

Quantum parameter estimation usually involves 
schemes for optimally estimating properties of quantum 
states or processes. The problem of parameter estimation 
for quantum processes can be treated as a task of opti¬ 
mal estimation of quantum states, where the states under 
investigation are the parameterized output states of the 
quantum processes. The estimation of unitary quantum 
processes is of interest and has been investigated in many 
different setups [25, '261. In a realistic scenario, the sys¬ 
tem evolution is unavoidably accompanied by some noise 
and the quantum process is treated as a general trace¬ 
preserving completely positive map, which is also called 
a quantum channel. The optimal estimation of a gen¬ 
eral one-parameter channel has been studied and some 
optimal estimation schemes have been proposed iza Eg. 
These schemes mainly consist of the identification of an 
optimal input state and an optimal estimator of the out¬ 
put state. In this paper, we restrict our investigation to 
unitary quantum processes, for simplicity, and concen¬ 
trate only on finding the optimal input states and the 
corresponding accuracy limits. 
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It is well known that the inverse of Fisher information 
gives the lower bound of the accuracy limit. Fisher in¬ 
formation characterizes the amount of information about 
the true value of a parameter that can be extracted from 
a probability distribution. Quantum Fisher information 
(QFI), which is a key concept in quantum metrology, is 
defined by maximizing the Fisher information over all 
possible measurements e og. The quantum Cramer- 
Rao theorem asserts that the precision is bounded from 
below by the inverse of the QFI Enusu- Due to its great 
importance in quantum metrology and parameter esti¬ 
mation, the QFI has attracted a lot of attentions EM3- 
However, analytical expressions for the QFI are difficult 
to obtain in most cases. 

The QFI associated with a state p for a parameter 0 is 
defined as [291 RT] 

F = Tr(pL 2 ), (1) 

where Tr stands for trace and L is the symmetric logarith¬ 
mic derivative operator, which is determined by dgp = 
(Lp + pL) /2. Consider a general unitary parametriza¬ 
tion process U = exp (—itH) for an initial pure state | ip), 
with the time-independent Hamiltonian H depending on 
a parameter 6, the QFI is simply connected with the vari¬ 
ance of a Hermitian operator PL in \ip), that is [42j [43] , 

F = 4((H 2 )-(n) 2 ), (2) 

where 

PL = i ( d e U t) U (3) 

is a Hermitian operator associated with the parametriza¬ 
tion. 

It has been shown that the variance of an Hermitian 
operator PL is maximized when the initial state | ip) is an 
equally weighted superposition of the eigenvectors |A m) 
and |A m ), which correspond to the maximum and mini¬ 
mum eigenvalues A m and A m of the operator PL, respec- 
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tively [U 02]; i-e., 

|^) = -^(|A M )+e^|A m )), (4) 

where (f> is an arbitrary relative phase, and the maximal 

QFI (MQFI) is 02 ] 

T max = (Am - A m ) 2 . (5) 

Here \ip) is the optimal input state. The convexity of 
the QFI guarantees that the MQFI is attained by a pure 
state rather than mixed states [44j. An important obser¬ 
vation is that \ip) can be a maximal entangled state for 
a many-body system, e.g., GHZ states, or NOON states. 
This is consistent with previous research where maximal 
entangled states are optimal input states to enhance the 
estimation precision. 

Previous research has mainly focused on the estima¬ 
tion of an overall multiplicative factor of a Hamilto¬ 
nian H (7j 0 HOj [El mi [38l Ho]. In such cases, the % 
operator is easily obtained as 'LL = —tdgH [03] ■ In a 
recent work |42| , Pang et al. studied a general Hamil¬ 
tonian parameter estimation problem where the param¬ 
eter does not appear as an overall multiplicative factor. 
Their study extends the quantum metrology to a more 
general case. Based on their work, we offered another ap¬ 
proach to get the QFI for unitary parametrization pro¬ 
cesses in [03]. In this approach, the QFI is related to 
the spectral decomposition of the initial state and the T-L 
operator of the parametrization process 

0 ° / - i\k- 1-1 

U = -td 9 H + iJ2±j->^H* k (deH), ( 6 ) 

where the superoperator H xk denotes a fcth-order nested 
commutator operation, H xk (-) = [H, ■ ■ •, [H 1 •]]. When 6 
is a multiplicative factor, as in the Mach-Zehnder inter¬ 
ferometer setup, H = 6 k H',{k = 1, 2, ...)(T2] 03], where 
H' is independent of the parameter to be estimated, the 
summation vanishes, and only the first term on the right- 
hand side of Eq. © contributes. 

In this paper, we investigate a general unitary 
parametrization process governed by an su(2) Hamilto¬ 
nian in the form of r • J, where r is a time-independent 
function of the parameter and J is the generator of su(2) 
algebra. This Hamiltonian covers many interesting ap¬ 
plications in physics. We derive the analytical expres¬ 
sion of the MQFI and find that it can be divided into 
two parts; one is quadratic in time t due to the depen¬ 
dence of the norm |r| on the parameter 0, and the other 
oscillates with time t because of the dependence of the 
direction e r = r/|r| on the parameter 9. Furthermore, we 
apply the theoretical expression in three practical exam¬ 
ples. The first two examples undergo time-independent 
su(2) parametrization, and the third one is subjected to 
time-dependent su(2) parametrization. 

This paper is structured as follows. In Sec II, we de¬ 
rive the H operator and the MQFI for a general su(2) 


Hamiltonian. In Sec III, three applications of the the¬ 
oretical result are listed. The first two are governed by 
time-independent su(2) Hamiltonians. We discuss the 
corresponding MQFI for the estimated parameters. In 
the third application, we further study the parameter es¬ 
timation with a time-dependent su(2) Hamiltonian and 
give the expression of the T~L operator and the MQFI. 
Finally, we give the discussion and conclusion in Sec IV. 


II. MAXIMAL QUANTUM FISHER 
INFORMATION FOR A GENERAL 
TIME-INDEPENDENT SU(2) HAMILTONIAN 

In this section, we investigate the 7~L operator and the 
MQFI for a general su(2) parametrization process. As¬ 
sume that the Hamiltonian takes the form of 

H = r • J, (7) 

where J = (J x , J y , J z ) is the generator of su(2) algebra 
and r = r (9) is a time-independent parametric curve in 
the parameter space. We denote the derivative of H with 
respect to 9 as 


d e H = v■ J, (8) 

where we define the velocity v = dr/d9 , which denotes 
the change of r over the parameter space of 9 , including 
both of the change of |r| over 9 and the change of the 
unit vector e r = r/|r| over 9. Thus the velocity v can be 
decomposed as v = v r + v t , with the radial velocity v r 
and the transverse velocity v t reading as 

v r = ( 9a ) 


v t = r 


, de r 
~d9 


(9b) 


Utilizing the commutation relation for su(2) algebra 
[a ■ J, b • J] = i( a x b) • J, the first two commutators in 
Eq. ([6| read 


H x dgH = i(r x v) • J, 
H x2 dgH = — [r x (r x v)[ • J 

= [I r | 2 v — (r • v)r] • J, 


(10a) 

(10b) 


where the relation a x (b x c) = (a • c)b — (a • b)c is 
used in the last equality. It is shown that both H x dgH 
and H x2 dgH are eigen operators of the superoperator 
H x2k , i.e., 


H x(2k+1) dgH = \r\ 2k H x d e H, (11a) 


H x{ 2 k+ 2 ) d g H = |r| 2 k H x 2 dgH, 


(lib) 
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where k = 1,2, .... By plugging Eq. (111 into Eq. ((6]), we 
have 


% = —tdgH + i ^ 


{it) 


2fc+2 


k—0 


(2k + 2)! 


rx(2 k +l)( de H) 




k =0 


{2k + 3)! 


= —tdgH + i (cos (| r1 1) — 1) 


H x dgH 


— (sin (|r|t) - |r|t) 


H x2 d g H 


( 12 ) 


and inserting the first two commutators in Eq. (101 into 
Eq. (12), we find that 

n = (r-v) (sin (|r|f) — |r| t) ^ ^ _ sin (|r|t) ^ ^ 


1 — cos (|r|t) 


(r x v) • J. 


(13) 


According to Eq. (|9|, we have that (r ■ v)r = |r| 2 v r , 
r x v = |r|e r x v t = |r|(v r x v t )/(d|r|/ci0); therefore, 
H can be further written as 


n = a j 

1 — cos(|r|<) 

\ r \M 

I 1 ! d8 


, . , , sin(lrlt) 

(V r X V { ) • J - fV r • J - |^|— -v t ■ J 

(14) 


where 


1 — cos(lrlf) . . sin(lrlt) . . 

A = I ' K x V ‘) - (15) 

l r l d9 1 1 


The three terms in Eq. (15) are perpendicular 
to each other and the norm of A is |A| = 

\Jv 2 t 2 + 4 ^ 5 -sin 2 therefore according to Eq. (JsJ) , 

the MQFI over the parameter 9 is 


F max = [j|A| - ( j)|A|] 2 


= 4 r 


v 2 t 2 + 4-^- sin 2 

r 2 


(16) 


where j is the maximal eigenvalue of J z . Equation (16) 
shows that the MQFI can be divided into two parts. The 
first part is quadratic in time t and is proportional to 
the square of the radial velocity v r . The second part 
oscillates with time t and is proportional to the ratio of 
the square of the transverse velocity v t to the square of 
r. According to Eq. (9b I, this ratio is equal to the square 
of the derivative de r /dd. Thus the quadratic term is due 
to the dependence of the norm |r| on the parameter 9, 
while the oscillation term is due to the dependence of the 
direction e r on 9. 

When t is small, we can expand Eq. (16) in t to the 

(17) 


second order, and the MQFI is simplified to 


F max = 4 j 2 t 2 [vj! + V* ] = 4 j 2 t 


2,_2 


That is, the MQFI grows quadratically with time t and 
depends only on the norm of v. 

If the change of |r| over the parameter space is 0, i.e., 
v r = ^e r = 0, then the MQFI is reduced to 

F max = 16 j 2 ^J sin 2 . (18) 

In this case, the MQFI always oscillates with the time 
so that its value is bounded. At the time points t = 
{2k + l)7r/|r|, k = 0,1,2,..., the MQFI reaches the opti¬ 
mal value, 


pmax 

-^op 



(19) 


III. APPLICATIONS 

In this section, we apply the theoretical expression 
in three practical examples. The dynamics of the first 
two examples are governed by time-independent su(2) 
Hamiltonians, and the third one is governed by a time- 
dependent su(2) Hamiltonian. Also note that the Hamil¬ 
tonians may depend on several parameters, and we only 
consider the single parameter estimation in this paper, 
i.e., we suppose that other parameters are known when 
we estimate a specific parameter. 


A. Time-independent Hamiltonians 


Case 1. Let us consider the Hamiltonian, ([7|, with 
in the explicit form in the spherical coordinates, 


r = r (sin 9 cos ip, sin 9 sin ip, cos 9 ), 


( 20 ) 


where r represents the amplitude of the external mag¬ 
netic field and 9, ip denotes the direction of the 
field. Suppose 9 is the parameter we want to esti¬ 
mate. The velocity vector is readily obtained as v = 
r (cos 9 cos <p, cos 9 sin ip, — sin 9). It is obvious that v r = 
0. According to Eq. (181, the MQFI over 9 is 


Fr*= 16 j z sin 2 (| 


( 21 ) 


and when t = {2k + l)n/r, k = 0,1, 2,..., F , ™ ax reaches 
its optimal value of 16j 2 . 

Similarly, if tp is the parameter to be estimated, the 
MQFI can be derived as 

E” ax = 16j 2 sin 2 0sin 2 (^^, (22) 

and when t = {2k + l)ir/r,k = 0,1,2,..., and 9 = 
pma.x reac } ies jt s optimal value of 16j 2 . In both cases, 
the MQFI is bounded. 
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FIG. 1. (Color online) The MQFI F^/ x as a function of time 
t with A = 1 and (a) wo = 0.1,(b) wo = l,(c) u>o = 10. 
Dashed (red) lines denote the contribution of the first term 
in Eq. (25 ) and solid (blue) lines denote the contribution of 
both terms, (d) We set wo = 1 and t = n/K , and F^ ax 
decreases as A grows. Here we take j = 1. 


Finally, if r is the parameter to be estimated, then 


= v 2 = 1 and the second term in Eq. (161 vanishes. 


Thus the MQFI reads 


F r max = 4j 2 t 2 . 


(23) 


In this case, the MQFI is unbounded and grows quadrat- 
ically with time t. 


Case 2. Let us consider an ensemble of N two-level 
atoms interacting with a static magnetic field. The 
Hamiltonian of this system can be written as 

H = ojqJ z + XJxi (24) 


where ujq is the atomic transition frequency and A is the 
Rabi frequency, which is proportional to the amplitude 
of the static magnetic field. 

Here we find the MQFI of ojq first. Comparing with the 
previous section, the coefficient vector is r = (A, 0 ,Wo) T . 
Denote the norm of |r| as K = yj A 2 + Wq; then the 
velocity is v = dr/duo = (0,0,1) T , the radial veloc¬ 
ity v 2 = d\r\/du)Q = ujq/K 2 , and the transverse veloc¬ 
ity v 2 = v 2 —v 2 = A 2 /K 2 , thus the corresponding MQFI 
is readily obtained as 


Fmax _ 4 -2 

uo J 


^0 4.2 

K 2 




(25) 


When t is very large, the oscillating term can be ne¬ 
glected. However, this term may also be important in 


some cases. We plot the MQFI .F“ ax with and without 
the oscillating term in Fig. [l] In Fig. [lja) , where we take 
A = 10wo = 1, we see that the effect of the oscillating 
term is very strong. In Fig. [ljb), we take A = Wo = 1, 
and the amplitude of the oscillating term is extremely 
small in comparison with the quadratic term. And in 
Fig. 0 c), where we take A = O.lwo = 1, the oscillating 
term could be neglected even when t is very small. 

On the other hand, in Fig. 0d), we set wo = 1 and t = 
n/K, and we can see that F 1 ™ 8 * decreases with increasing 
A. This can be explained as follows. When A is very large, 
the contribution of the first term ujqJ z in Eq. (24) can be 
neglected; i.e., when A wo, i^ ax —> 0. 

Following a similar procedure, the MQFI with respect 
to A is obtained as 


f max _ 4^-2 


A 2 n 4 wq 
—prt 2 + — J- sin 
K 2 K 4 


2(?) 


(26) 


This result is similar to Eq. (251 and we omit the discus¬ 
sion here. 


B. Time-dependent Hamiltonian 

Hitherto all the Hamiltonians considered in this paper 
are time independent. If the Hamiltonian is time depen¬ 
dent, generally the unitary operator will involve a time¬ 
ordering procedure and Eq. © fails. In this subsection, 
we go beyond the previous situations and investigate a 
time-dependent Hamiltonian. In the following we replace 
the static magnetic field in the preceding subsection with 
a time-dependent driving field, and the Hamiltonian is 

H = oj 0 J z + A [J x cos(wf) + J v sin(wf)] . (27) 

Here w is the frequency of the driving field. This 
Hamiltonian could reflect many realistic physical sys¬ 
tems, such as in the Ramsey spectroscopy and NMR 
techniques era. The FL operator of this Hamiltonian 
cannot be written in the form of Eq. (J6]) and the previ¬ 
ous method must be changed. We can circumvent this 
problem by utilizing a rotating transform. In a rotating 
frame, the original state vector \i/>) is transformed as 

\$)=R\ 1 >) = e™ J »\il;), ( 28 ) 

and the effective Hamiltonian in this new frame is 

+ dB) 

H eS = RHrf - idi¬ 
ot 

= A J z + A J x , (29) 

which is time independent. Here A = wo — w denotes the 
detuning between the atom transition and the driving 
magnetic field. 

Therefore, the unitary evolution in the original frame 
is readily obtained as 

U = U 1 U 2 = e ~ iutJz e~ iHef[t , 


(30) 
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with U-[ = W = e~ lultJz and U 2 = e~' lHeSt . In such away, 
we turn the time-dependent Hamiltonian of Eq. (27) into 
two time-independent ones, Hi = u>J z , and The H 
operator for two consecutive unitary operation U = U\ U 2 
can be derived as 

H = ide C/1C/2 

= H 2 + u\UiU2, (31) 

with Hi = idgU\Ui and H 2 = id gUl Uz . Here 0 can be 
A,wo, or uj. For the case in Eq. (301, if A or loq is the 
parameter to be estimated, it is obvious that H 1 = 0 and 
one can simply get the MQFI by replacing wo in Eq. (25) 
and (26) with A . When A = 0, F” 
value of 4j 2 f 2 . 


71 max 
A 


attains its optimal 


Next, let us find the H operator with respect to oj. 
The Hamiltonian associated with H 2 is H 2 = H e ff, and 
the corresponding coefficient vector reads r = (A, 0, A) T . 
Then by utilizing Eq. © , we can obtain that 


W = ^72 (-A,0,-A), 
—AA , 

Vr X Vi = (0, 1,0) . 


(32) 


(33) 


K' 2 

Therefore, according to Eq. ( fl4| , we have that 

H 2 = -^3 [AA [K't - sin (K't)] J x 

+XK' [1 - cos(A"Q] J y 
+ [A 2 K't + A 2 sin {K't)] J z ]. 

On the other hand, the Hamiltonian associated with 
Hi reads Hi = coJ z , where the parameter enters in the 
Hamiltonian as a multiplicative factor. According to 
Eq. © , we have Hi = — tJ z . Based on the formula 

exp(A)B exp(—A) = exp(A x )F = X)fclo NTT- 8 , we ob “ 
tain 

U\HiU 2 

= — [AA [1 — cos(A'f)] J x + \K' sin (K't)J v 

+ [A 2 + A 2 cos {K't)] J z \ , (34) 

where we denote K' = \/A 2 + A 2 . 

Combining Eqs. (31), (331, and (341, we have 

= lsill(A " i) HAAA3 { _ Aaj , + A v.) 

[1 — cos (K't) — K't sin (K't)] 


K' 2 


-XJy. 


(35) 


Based on the norm of the vector A', the MQFI is readily 
obtained as 


A 2 


= 4/^71 [2 + K ' 2 t 2 

- 2 K't sin (K't) - 2 cos (K't)] . (36) 



A At 


FIG. 2. The MQFI F’™ ax as a function of the detuning A 
(a) and A (b). (a) We take A = 1, and t = 1. F’™ ax attains 
the optimum when A = 0. (b), We take A = 0 and A = 1. 
F™ ax grows quadratically with time t. We set j = 1 in both 
panels. 


Figure 2.(a) shows that when A = 0 (resonance), F™ ax 
attains its optimal value. This can be proved theoreti¬ 
cally by taking the derivative of Eq. ( |36| , 

^TT ax |A=o=0. (37) 

When A = 0, we plot _F™ ax in Fig. 2(b) with A = 1. As 
t increases, the MQFI can be approximated as _F™ ax = 
4j 2 t 2 and the oscillation can also be neglected. 


IV. DISCUSSION AND CONCLUSION 

In the previous sections, we have only discussed the 
single-parameter estimation problem. The extension to 
multi-parameter estimation is straightforward and has 
been discussed in our recent work | l%5] . 

The quantum Cramer-Rao theorem asserts that the 
variance of any unbiased estimator is bounded below 
by the inverse of the QFI. In this paper, we compute 
the MQFI for su(2) parametrization processes, thus get¬ 
ting the corresponding ultimate estimation bound. The 
achievement of the bound is of practical interest and var¬ 
ious methods of building the optimal estimators are dis¬ 
cussed in the literature, including Bayes estimators and 
maximum-likelihood estimators 01 HZj. The optimal 
quantum estimator is also discussed in |29| and 05| for a 
general process. The specific form of the optimal estima¬ 
tor for an su(2) parametrization process will be discussed 
in our future works. 

In summary, we investigate a unitary parametrization 
process governed by an su(2) Hamiltonian r • J. We find 
the optimal input state and the MQFI when the Hamil¬ 
tonian is independent of time. The optimal input state 
\ijj) in Eq. © can be a maximal entangled state for a 
many-body system. A similar result is discussed in []0], 
where the author found that the optimal input state for 
an isotropic depolarization channel is also a maximally 
entangled state. We show that the MQFI can be divided 
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into two parts. One part is quadratic in time t due to 
the dependence of the norm |r| on the parameter 9, and 
the other part oscillates with time t because of the de¬ 
pendence of the direction e r = r/|r| on the parameter 9. 

We apply this result to a typical scenario and find the 
MQFI corresponding to the amplitude and the direction 
of the magnetic field. We also investigate the MQFI for 
an ensemble of atoms interacting with a static field. We 
find that in some cases the oscillating terms can con¬ 
tribute significantly and thus can not be neglected. We 
further investigate a time-dependent driving field and 
find that the MQFI of the field frequency reaches its op¬ 
timum when the driving field is in resonance with the 


atomic transition frequency. Moreover, the MQFI of the 
amplitude of the external field also reaches its optimum 
when the field is in resonance with the atoms. 
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